■ APPROXIMATIONS OF THE SPECTRAL FUNCTION* 

O 



o 



< 



oo 
o 

O 



ARTUR M. ANKOWSKI 



|5Pj Institute of Theoretical Physics, University of Wroclaw, 

pi. M. Borna 9, 50-204 Wroclaw, Poland 



artank@if t . uni . wroc . pi 



>. 

The ICARUS and future liquid argon neutrino experiments generate 
demand for evaluating the spectral function of argon. In this paper we use 
oxygen nucleus as a testing ground for our phenomenological approach to 
the spectral function and probe the influence of momentum distribution 
and treatment of the mean field spectral function on the differential cross 
sections. The obtained model reproduces very well results of the exact 
spectral function of oxygen and can be applied to heavier nuclei, such as 
calcium or argon. 



PACS numbers: 13.15.+g, 25.30.Pt 
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1. Motivation and outline of the paper 

^ ' Thanks to experience already gained with the ICARUS T600 TPC [Q 

one knows that liquid argon (LAr) has many advantages in neutrino exper- 
iments. They make LAr technology interesting for planned detectors, e.g. 
for T2K [2] and NuMI 3 . To use them fully, one has to reduce uncertainties 
by evaluating nuclear effects as precise as possible. From (e, e') scattering 
it's clear that the Fermi gas model can only be the first approximation. 

More elaborate approach is the spectral function (SF) formalism. Where 
the impulse approximation is valid (i.e. when neutrino energy E u is greater 
then a few hundred of MeV §1]), the SF describes nuclei most accurately. 
It was applied to (I, I') scattering previously [3]- The problem is that the 
exact SFs exist only for a few double closed shell nuclei, 3 He, and nuclear 
matter (HI HE]. For argon fgAr exact computations cannot be performed, 
so one is forced to seek for the best possible approximation. Opportunity 
to verify the quality of given approximation is provided by oxygen nucleus, 
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where the exact SF exist {2]: applying the same method to calculation of 
the oxygen SF and comparing the result to the exact SF can give a notion 
of discrepancies between them. 

The first attempt at constructing the argon SF is Ref. [S]. Presented 
there approximation gives the differential cross sections which somewhat 
differ from the corresponding ones for the exact SF (compare Fig. In 
the paper presumption was made that the discrepancies come from: 

- oversimplified treatment of the mean field spectral function, 

- different momentum distributions in the spectral functions. 

This work is devoted to detailed study of the two effects. After a brief 
introduction into the SF approach in Sec. [2] and describing the simplest 
approximation of the SF in Sec. 01 (i.e. the approach of [§]), we consider in 
Sec. H] an influence of NN-correlations on the mean field SF. Then in Sec. El 
the sensitivity of the SF on the momentum distribution is discussed. 

Conclusions allowed for working out in Sec. El a satisfactory approxima- 
tion of the exact spectral function of oxygen, which can be applied to argon 
and other nuclei. 



2. Basic information on the spectral function 

The spectral function (SF) of a given nucleus P(p, E) is the probability 
distribution of finding a nucleon with momentum p and removal energy E. 
Formal definition can be expressed as jlUj 

P(p,E) = (i(M A )\a\p)5(H - M A + M - E)a(p)\i(M A )} , (1) 

where H is the intrinsic Hamiltonian of the (A — l)-nucleon system, \i(M A )) 
denotes the state of the initial nucleus of mass M A (assumed to be at rest), 
and M is nucleon mass. 

Oxygen nucleus is isoscalar (N = Z), therefore its SF consists of only 
two parts: the mean field and the short-range correlation one: 

P(p, E) = [P MF (P, E) + P corr (p, E)] . (2) 

All the knowledge about low-energy and low-momentum nucleons (i.e. 
the shell model information) is included in the mean field term Pmf(Pj-E')- 
Corrections to the independent-particle behavior of nucleons are described 
by the correlation part P C orr(P; E). From the nuclear matter calculations one 
knows that at high energy and momentum such correlations are dominated 
by the two- nucleon interactions [TT] . 
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The mean field SF consists of contributions of every shell model state a 
below the Fermi level of |12j : 

P MF (p,E) = j Yl c a A a \^ a {p)\ 2 5(E a + E R { v )-E). (3) 

In the above equation E Q is the energy of the state a described by the single- 
nucleon wave function 4> a (p) (normalized to 1) with the occupation proba- 
bility c a , and the number of particles A a , whereas Er(p) = p 2 /(2M^i) is 
the recoil energy of the residual nucleus. 

Interactions between nucleons cause their partial redistribution from the 
states below the Fermi level to the levels of higher energy, and occurrence 
of non-zero width of each level. It means that when NN-correlations in the 
MF part are taken into account, delta function in Eq. (J3J) should be replaced 
by appropriate distribution |7j. We shall return to this issue in Sec. EJ 

Due to Eq. the momentum distribution of nucleons reads 

n(p) = {i{M A )\a\p)a{p)\i{M A )) = j P(p,E)dE (4) 
and as a consequence of Eq. © consists of two contributions ^2] : 

n (p) = 1~ Kmf(p) + ^corr(p)] • (5) 



3. Simplest approximation of the SF 

Firstly note that Eq. (j3J) simplifies significantly when one replaces E a 
in the argument of the delta function by the average separation energy 

E {1) =E a<aF Co l E a A a /Am-- 

P$W(p,E) = n MF (p)5{EW +E R (p) -E), (6) 
because then only the momentum distribution 

«mf(p) = Yl c*\Mp)\ 2 A a /A (7) 

a<ap 

occurs, which is averaged over single-particle levels. 

Secondly, as it was mentioned in the previous section, major contribution 
to the correlations between nucleons comes from the two-nucleon process. It 
consists in forming a cluster by two nucleons with high relative momentum 
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Fig. 1. Differential cross section da/dE^ of quasielastic scattering off O ob- 
tained from the Fermi gas model (dotted line), the Benhar's spectral function 
(dashed line) and the simplest approximation of the spectral function 9 (solid line) . 
Left: Pauli blocking is absent. Right: Pauli blocking included. 



while the other (^4 — 2) nucleons remain soft |13| . When we restrict ourself 
to such an interactions, the correlation SF can be expressed as |14j : 

P corr (p, E) = n coxr (p) 1 — j-4 /- [exp(-ap min ) - exp(-ap max )] . (8) 
iPl V 7r 

Occurring here a = 3/(4(p 2 )/3) is inversely proportional to the mean value 
of the MF momentum squared (p 2 ) times = (A — 2) /(A — 1) and 

0) 



2 

Pmin 



0\p\ - V2M[3[E - W> - E R {p)} 



P 2 



/3|p| + V2M(3[E - W> - E R { V )) 



The threshold value E^ is interpreted as the two-nucleon separation energy 
averaged over low-energy configurations of the (A — 2)-nucleon system and 
can be approximated by E^ = Ma-2 + 2M — Ma- 

We will refer to the described approach as "the simplest approximation 
of the spectral function" or "the simplest spectral function" (SSF). 

Figure ^ illustrates that the cross section dcr/dE^ of both SFs clearly 
differs from the one of the Fermi gas. As long as Pauli blocking (PB) is not 
included, the SSF reproduces the result of the Benhar's (i.e. exact) SF quite 
well. Presence of PB enlarges difference between them, because it changes 
the shape of the peak in other way. 

In the next section we investigate whether more elaborate approximation 
of the mean field SF removes this discrepancy. 



Ankowski printed on February 9, 2008 



■5 




E (MeV) 

Fig. 2. Our fit of the Gaussian distribution of the level \s\ (solid line) to appropri- 
ate region of the energy distribution Smf(E) calculated from the Benhar's spectral 
function of 16 O (dotted line). 



4. Improved treatment of the mean field SF 

In Sec. |2]we mentioned that NN-correlations broaden energy levels, so in 
Eq. J2J) instead of delta function an appropriate distribution should occur [7j : 

P MF (P,£) = jJ2 c a A a \(t> a {p)\ 2 F a {E a + E R {p) - E). (10) 

Q<ap 

In order to simplify it, let's replace the single-particle momentum distribu- 
tion |0 Q (p)| 2 by the average riMF(p) defined in Eq. (0: 

(p)-r Yl c a A a F a (E a + E R (p) - E). (11) 
We decided to use the Gaussian distribution 



F a (x) = l/(V^D a ) exp [- {x/D a y\ . (12) 

To describe a level a we need to now its difiuseness D a and its "mean 
field" occupation probability c a . If one used the "raw" c Q 's (given for ex- 
ample in J3])) partial double counting would occur due to contribution of 
the correlation term of the SF (raw value says only that the nucleon can be 
found in given state with certain probability, not how it get there: whether 
it happened due to "natural placement" or due to "correlation kick" ) . 

For 16 we shall determine the value of D a from the Benhar's SF, for 
other nuclei one should get them from direct calculations (e.g. for loCa the 
data is available in JB])- Oxygen is a quite heavy nucleus in the sense that 
it's recoil energy for the mean value of momentum |p| = 180 MeV/c [H] is 
~ 1 MeV. When we neglect it, the separation 



P M f{p,E) = n M F{p)SMF{E). 



(13) 
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Fig. 3. Sensitivity of the quasielastic differential cross section da/dE^ on treatment 
of the mean field spectral function. The simplest (dotted line) and the Gaussian ap- 
proximation (solid line) compared to the Benhar's spectral function. Left: Results 
without Pauli blocking. Right: Pauli blocking included. 



holds, and we easily get that 

S M f(E) oc J P MF (P, E) d'p = J [P(p, E) - P corr (p, E)) d 3 p. (14) 

We inserted in the above relation the Benhar's SF and calculated the en- 
ergy distribution Smf(E). By fitting Eq. l|12j) successively to every peak 
corresponding to the energy level, we determined all D^s, see Figure HI 

Note that Smf(E) shown in Fig. [2] vanishes for E > 87 MeV. For oxy- 
gen P C orr(p ; E) appears above the threshold energy = 26.33 MeV 
higher then the energy of quite sharp levels lp| and lp^, hence the cor- 
relations admixes only to c a of ls^ level. After subtracting calculated value 
of J S COII (E)dE over E € [26.33; 87] MeV from c a of lsi, we get the proper 
"mean field" value. Now we have all the indispensable parameters to obtain 
predictions for the model. 

The cross sections for this approximation are shown in Fig. |31 the mo- 
mentum distribution is the same as in the simplest SF and so is the cor- 
relation part of the SF, the only change is the different treatment of the 
mean field SF. This approach causes that the shape of the cross section 
is similar to one of the exact SF, regardless of including PB or not. The 
height of the peak is also slightly reduced, therefore the agreement between 
the approximation and the Benhar's SF is now better. 



5. Dependence on the momentum distribution 

To study how the difference between the momentum distributions affects 
difference between the cross sections, we divided the Benhar's distribution 
into MF and correlated part as in Fig. |1J we kept the same value of n corr (p) 
for p = and took care of its smooth transition into the original distribution 
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Fig. 4. Contributions to the momentum distribution of nucleons: the mean field 
(dashed) and correlation part (dotted line) sum up to the total momentum distribu- 
tion (solid line). Left: Distribution from |12j. Right: Corresponding one calculated 
from the Benhar's SF, by analogy divided into two parts (details in text). 
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Fig. 5. Sensitivity of the quasielastic differential cross section do~/dE n on applied 
momentum distribution. The simplest approximation with the momentum distri- 
bution used previously 9 (from ^21; dotted line) and with "the Benhar's distri- 
bution" (see Fig. 0J solid line) compared to the Benhar's spectral function (dashed 
line). Left: Without Pauli blocking. Right: With Pauli blocking included. 



for high values of p. The calculated components of the total distribution 
were applied to the simplest approximation of the spectral function (SSF). 

As it follows from Fig. nearly whole discrepancy between the cross 
section for SSF obtained in this way and the corresponding one for the 
Benhar's SF disappeared when Pauli blocking (PB) is absent. Inclusion of 
PB reveals different behavior of the two approaches. 

We conclude that the momentum distribution is responsible for the 
height of the peak, but not for it's shape (compare solid and dotted line in 
Fig. 03). Therefore the SSF diverges from the exact SF when PB is present. 

Another important information is the high sensitivity of da/dE^ on 
n(p). This feature is quite inconvenient, because it suggests that one should 
take into account the difference between momentum distribution of protons 
and neutrons, if appropriate prediction of da/dE^ for high En is required. 
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Fig. 6. Gaussian approximation with "the Benhar's distribution" (solid line) is 
significantly better then the simplest approximation with the original distribution 
(dotted line) and nicely reproduces results of the Benhar's spectral function (dashed 
line). Left: Without Pauli blocking. Right: With Pauli blocking. 

6. Combination of the two effects 

As we have seen in Sec. |1J the enhanced treatment of the mean field 
SF affected the shape of the cross section's peak, but it was a bit too high, 
whereas the change of the momentum distribution in Sec. |S] reduced the 
hump's height, however it's shape disagreed. When we combine the two 
improvements, using both "the Benhar's momentum distribution" and the 
Gaussian mean field SF, discrepancy vanishes almost completely, see Fig- EI 

We conclude that the remarks made in [3] were right: the cross section 
of the Benhar's SF can be satisfactory reproduced when one uses the same 
momentum distribution and the refined treatment of the mean field SF. 

Our approximation remains simple enough and can be applied to other 
nuclei. The next goal will be loCa where the needed nuclear data is available. 

The author would like to express his gratitude to Jan T. Sobczyk for 
stimulating discussions on the spectral function and to Omar Benhar for 
providing his spectral function of oxygen. 
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